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ABSTRACT 

The final mass of a newborn star is set at the epoch when the mass accretion onto the star is terminated. 
We study the evolution of accreting protostars and the limits of accretion in low metallicity environments. 
Accretion rates onto protostars are estimated via the temperature evolution of prestellar cores with different 
metallicities. The derived rates increase with decreasing metallicity, from M ~ 10"^ Mq/yy at Z = Zq to 
10"^ Mq/yt at Z = 0. With the derived accretion rates, the protostellar evolution is numerically calculated. 
We find that, at lower metallicity, the protostar has a larger radius and reaches the zero-age main-sequence 
(ZAMS) at higher stellar mass. Using this protostellar evolution, we evaluate the upper stellar mass limit 
where the mass accretion is hindered by radiative feedback. We consider the effects of radiation pressure 
exerted on the accreting envelope, and expansion of the H II region. The mass accretion is finally terminated 
by radiation pressure on dust grains in the envelope for Z > 10"^ Zq and by the expanding H II region for lower 
metallicty. The mass limit from these effects increases with decreasing metallicity from M* ~ 10 Mq at Z = Z© 
to ~ 300 Mq at Z = 10"^ Z©. The termination of accretion occurs after the central star arrives at the ZAMS at 
all metallicities, which allows us to neglect protostellar evolution effects in discussing the upper mass limit by 
stellar feedback. The fragmentation induced by line cooling in low-metallicity clouds yields prestellar cores 
with masses large enough that the final stellar mass is set by the feedback effects. 

Subject headings: accretion - stars: early-type - stars: evolution - stars: formation- stars: pre-main-sequence 



1. INTRODUCTION 

What determines the mass of a newborn star? Our knowl- 
edge on this fundamental question in astrophysics is still lim- 
ited. Observationally, the stellar initial mass function (IMF) 
is a power law with the Salpeter slope -2.3-2.6 on the high- 
mass side and becomes flatter below a break at ~ IM© on the 
lower-mass side. Towards even lower mass (< O.IM©), the 
number of stars appears to decline (e.g., Muench et al. 2002). 
That is, the characteristic mass in the solar neighborhood is 
~0.1-1Mq. 

Recently, the mass distribution of dense cores in a star- 
forming region was found to have a similar shape to the stellar 
IMF, with a mass offset about a factor of three (Alves et al. 
2007). Also, the typical mass of the cores Mcore is approxi- 
mately the Jeans mass Mj: 

/ « \-i/2/ T 
M..^M,.2M©(^^) (^^j , (1) 

where we used the typical number density n ~ 10^ cm""* and 
temperature of the cores T ~ 10 K. These facts indicate that 
the characteristic stellar mass is set akeady at the formation 
stage of dense cores through gravitational fragmentation. The 
factor of three reduction of the stellar mass from the core mass 
can be a result of some stellar feedback, e.g., protostellar out- 
flows, which reduce the amount of material in the core, or 
binary/multiple formation. Theoretical estimations support 
such an efficiency for the low-mass cores (e.g., Matzner & 
McKee 2000). 

Even though the characteristic stellar mass is determined by 
the gravitational (or turbulent) fragmentation process, forma- 
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tion of the most massive stars will be limited by the feedback 
from the protostar, as inferred from the following considera- 
tion. A small protostar, which forms at the center of the core 
subsequent to its gravitational collapse, grows in mass by ac- 
cretion from the leftover core material, i.e., envelope. Since a 
massive protostar becomes so luminous during accretion, ra- 
diation pressure on dust grains easily disrupts the accretion 
envelope (Larson & Stan-field 1971; Kahn 1974; Wolfire & 
Cassinelli 1987). Under a typical accretion rate, i.e., a Jeans 
mass Mj of gas falling on the protostar in a free-fall time fff, 

M^^^-.2xl0-M©/yr(^j ' 

where is the sound speed in the natal core and G is the grav- 
itational constant, the upper mass limit by radiative pressure 
falls on about 10-20M©. 

Although rare, stars more massive than this limit do exist in 
the Milky Way. How such stars formed is still debated. Ac- 
cretion through a circumstellar disk can somewhat alleviate 
the radiative feedback. Besides, it appears that massive star 
formation requires a non-standard environment. For exam- 
ple, faster accretion than the standard or stellar mergers are 
sometimes invoked. In fact, observations tell us that accretion 
rates on massive protostars are quite high, M > 10"^ Mg/yr 
(e.g., Zhang et al. 2005; Keto & Wood 2006; Grave & Ku- 
mar 2009), and the natal dense cores are unusually massive 
McoK > 10 M© 10 Mj (e.g., Saito et al. 2006; Motte et al. 
2007). However, as a guide, the upper mass limit within the 
standard accretion scenario is still important. 

Metallicity governs the thermal balance in prestellar cores, 
thereby affecting the core temperature and the protostellar ac- 
cretion rate through equation (|2|. The accretion rate in turn 
alters the protostellar evolution. Metals also act as catalysts 
in the CN cycle, modifying the stellar structure. In addition, 
lowering the dust opacity is expected to reduce the strength of 
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Fig. 1 . — Top panel: Temperature evolution of prestellar cores with differ- 
ent metallicities (taken from Omukai et al. 2005). Tlie Jeans mass is constant 
along the thin dotted lines, Mj = 10"-, 1 and 10-, and 10* Mq. In calculating 
Mj, the hydrogen is assumed to be fully atomic (molecular) for A/j > 10-Mq 
(Mj < 1M0, respectively). Bottom panel: The fiducial accretion rates de- 
rived from the thermal evolution for different metallicities (see text). In both 
panels, the solid and dashed curves alternately represent the cases with metal- 
licities Z=0, 10"*, 10"\ 10"*, 10"^ 10"^ 10"' andZg. 

the radiative feedback. As a first step, we here aim to calcu- 
late the evolution of accreting protostars and evaluate the mass 
limit where the accretion is prohibited by the stellar radiative 
feedback as a function of metallicity under the standard accre- 
tion rates and spherical symmetry. We consider the effects of 
radiation pressure exerted on the accretion envelope, as well 
as expansion of an H II region. 

The organization of this paper is as follows. In § |2l we 
calculate the evolution of protostars at various metallicities. 
Next in § [3] we examine the radiative feedback effects. We 
also consider the limit of mass accreted over the main se- 
quence lifetime. The final stellar mass is eventually derived 
as a function of the metallicity. In § |4] we also compare the 
obtained final stellar masses with the typical core masses at 
each metallicity. §|5]is devoted to conclusions. 

2. PROTOSTELLAR EVOLUTION WITH VARIOUS METALLICITIES 

2.1. Accretion Rates onto Protostars 

The accretion rate onto the protostar is determined by the 
prestellar temperature (eq. |2]l, which depends on metallic- 
ity. The temperature evolution during the prestellar collapse 
for various metallicities has been studied by Omukai et al. 
(2005) by using a one-zone model with detailed chemical and 
radiative processes. The prestellar core is assumed to undergo 
run-away collapse, i.e., its central flat part has a radius about 



the Jeans length and the central density increases at a free-fall 
timescale. Note that its radius and mass generally decreases 
with increasing density and thus most of the material origi- 
nally in the core is left behind after some moment. We adopt 
the result by Omukai et al. (2005) in calculating the accre- 
tion rate, and their Figure 1 is reproduced in the upper panel 
of Figure [T] In general, prestellar temperatures are higher for 
lower metallicity. 

The thermal evolution track (p, T) is converted to the ra- 
dial density or temperature profile of the envelope (r = Aj,p 
or T), where Aj is the Jeans length. When materials within a 
radius r have accreted on the protostar, the protostellar mass 
is ~ Mj, where Mj is the Jeans mass with (p, T) at radius 
r. The accretion rate to this protostar is ~ c] /G, where 
is the sound speed at radius r. We derive accretion rates as 
a function of the instantaneous protostellar mass by convert- 
ing the thermal evolution track (p, T) to (M, = Mj,M = /G). 
The lower panel of Figure [T] shows those accretion rates. The 
higher accretion rates at lower metallicities reflect the higher 
temperatures in the envelopes. 

Here, we have adopted the standard accretion rate M ~ 
Cj /G. However, the accretion rate can vary according to the 
way of prestellar collapse. For example, in a dynamical run- 
away collapse, the Larson-Penston (LP) solution, the accre- 
tion rate reaches as high as M ~ 47 / G, while the stan- 
dard rate is almost exact for the static collapse of the singu- 
lar isothermal sphere, the Shu solution. In realistic hydro- 
dynamical calculations, the prestellar collapse asymptotically 
approaches the LP solution. The accretion rate is initially 
as high as ~ 10 times the standard rate and decreases grad- 
ually until it reaches the standard rate (e.g., Larson 2003). 
To accommodate possible such variations, we also examine 
the evolution with the accretion rate ten times higher than the 
standard value. 

2.2. Result of Protostellar Evolution 

Next, we calculate the protostellar evolution with the de- 
rived accretion rates at each metallicity. Stellar feedback ef- 
fects, which possibly alter the accretion rates, are separately 
examined in § 13.11 with the calculated evolution. Details of 
our numerical method are summarized in our previous pa- 
pers (e.g., Hosokawa & Omukai 2009 and references therein). 
We explain a brief outline of the method and some improve- 
ments. Figure |2] schematically depicts the structure of a pro- 
tostar and accretion envelope. Dust grains in the accretion 
envelope evaporate at a dust destruction front, where their 
temperature reaches about 2000 K. The envelope is divided 
into inner gas envelope and outer dust cocoon, whose bound- 
ary is the dust destruction front. We numerically solve the 
structure of an accreting protostar and gas envelope. The ba- 
sic equations for the protostar are the stellar structure equa- 
tions with mass accretion (e.g., Stahler, Shu & Taam 1980, 
Palla & Stahler 1991). We adopt the free-fall flow for the 
optically thin gas envelope. When the flow becomes opaque 
before reaching the stellar surface, we solve the structure of 
the opaque part, called a radiative precursor. The flow fi- 
nally hits the stellar surface and forms an accretion shock 
front. The whole structure of the protostar and gas envelope 
is consistently determined to satisfy the shock jump condi- 
tions. An evolutionary calculation begins with a small ini- 
tial protostar constructed according to Stahler, Shu & Taam 
(1980). We construct sequence of models with increasing stel- 
lar mass by accretion, until the protostar reaches the zero-age 
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Fig. 2. — Schematic view of a protostar and its surrounding accretion envelope. The accretion flow hits the stellar surface at the accretion shock front. If the 
flow becomes opaque before hitting the stellar surface, a photosphere emerges outside the stellar surface. The optically thick part of the flow is called a radiative 
precursor Dust grains in the envelope evaporate at a dust destruction front. The layer outside the destruction front is a dust cocoon. With the optically thick 
cocoon (left figure), most of the UV and optical light is absorbed by grains just outside the dust destruction front and is re-emitted as infrared light. With the thin 
cocoon (right figure), on the other hand, the stellar UV and optical light directly escapes out of the cocoon. 
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Fig. 3. — Evolution of protostellar radii with metallicities Z = 0, 10 , 
10^, 10"- and Zq. The upper and lower panels show the cases with the 
fiducial and ten times higher accretion rates, respectively. The open circles 
on the curves mark the epoch when the energy production rate by deuterium 
burning reaches 80 % of the steady burning rate, Ld st = M5d, where (5d is 
the energy available from the deuterium burning per unit gas mass. The filled 
circles on the curves indicate the epoch of ZAMS anival defined as when the 
total energy production rate by hydrogen burning reaches 80 % of the stellar 
luminosity L* . 



main-sequence (ZAMS) stage. Nuclear burning reactions of 
deuterium, hydrogen, and helium are included in our calcu- 
lations. Further evolution in the post main-sequence stage is 
out of scope of this paper. In § 13.21 we roughly estimate when 
the MS stage is over with mass accretion. In this paper, we 
use the tabulated equation of state generated following Pols et 
al. (1995). Opacity tables at < Z < Zq are also generated 
by logarithmic linear interpolation between the tables at Z = 
and Z = Zq . 

Figure [3] presents the calculated protostellar evolution at 
each metallicity. The evolution at Z = agrees with previ- 
ous calculations by Omukai & Palla (2001, 2003). There are 
four evolutionary stages with the fiducial accretion rates (up- 
per panel): (i) adiabatic accretion (M» < 6 Mq), (ii) swelling 
(6 Mq<M^< 10 Mq), (iii) Kelvin-Helmholtz (KH) contrac- 
tion (10 Mq < M„ < 100 Mq), and (iv) accretion to ZAMS 
star (M* > 100 Mq). A key quantity to understand varieties of 
the evolution is a ratio between two timescales; the accretion 
timescale, 

t =^ (3) 

'acc , • ; \-^/ 

M 



and the KH timescale, 



GMl 



(4) 



The two timescales face and Jkh can be regarded as evolution- 
ary and thermal-adjustment timescales of the protostar The 
protostar gains materials with high entropy generated at the 
accretion shock front on face- Meanwhile, radiation diffuses 
outward and reduce the stellar entropy on Jkh- In the early 
phase, face is much shorter than Tkh- The entropy taken into 
the star adiabatically accumulates in the stellar interior (adia- 
batic accretion). This stage is over because Tkh decreases (L* 
increases) as the stellar mass increases. As Jkh approaches 
face, the interior entropy is gradually transferred outward by 
radiation. A part of transferred entropy is received by the 
outer layer near the stellar surface. This causes the swelling 
of the star Finally, face becomes longer than Jkh- The star 
loses energy by radiation, and contracts hydrostatically (KH 
contraction). The central temperature increases with the KH 
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contraction. When hydrogen burning is ignited, the protostar 
reaches the ZAMS stage. The above four evolutionary stages 
can be appHed for other Z > cases. The protostellar radius 
gradually increases in the early stage and then experiences 
the swelling. The protostar contracts in the next stage and fi- 
nally arrives at the ZAMS. In all cases, the timescale balance 
switches from face < fKH to face > fKH at the swelling. However, 
Figure [3] also presents some different aspects of the evolution 
with different metallicities. The following three evolution- 
ary tendencies result from the lower accretion rates at higher 
metallicities (see also Hosokawa & Omukai 2009). 

(i) The protostellar radius is larger for lower metallicity. For 
Z = 0, the radius exceeds 100 Rq at Z = 0, while only a few Rq 
for Z = Zq. This reflects higher entropy content in the star for 
higher accretion rate (see e.g., Hosokawa & Omukai 2009) as 
the radius is related with entropy s by 



-1/3 

ocM, ' exp[consf. x s]. 



(5) 



Without enough time after the shock for accreted matter to 
cool radiatively before burying into opaque interior, the inte- 
rior entropy is higher for higher accretion rate, (ii) The pro- 
tostellar mass is higher at the onset of ZAMS phase for lower 
metallicity. The star arrives on the ZAMS at ~ 100 M© for 
Z = 0, and at a few Mq for Z = Z© (also see Fig|4]l. Since the 
stellar radius is larger for lower metallicity (e.g., lower accre- 
tion rate), the interior temperature is higher according to the 
relation, 

T = —— (X — . (6) 
Tip R^ 

Thus, higher is needed for ignition of hydrogen burning at 
~ 10^ K. (iii) Signatures of deuterium burning appear only at 
metallicity Z > lO'^Z©. At Z = Z©, the protostar becomes al- 
most fully convective by deuterium burning at M, ~ 0.1 Mq. 
The linear increase of the radius at 0.1 M© < < 0.3 M© 
is due to this effect. After that, the star gradually returns to be 
radiative by the swelling at M^, ~ 2 Mq. This is in contrast to 
the metal-poor cases. At Z = 0, deuterium burning begins later 
atM* ~ 10 M© and hardly affects the evolution. The protostar 
remains radiative until its arrival to the ZAMS. These varia- 
tions are also due to different evolution of the stellar inner 
temperature following equation (|6]l. 

Apart from different accretion rates, metallicity still influ- 
ences the protostellar evolution. For example, the accreting 
MS star has smaller radius at the lower metallicity. This is 
due to different initial abundances of C and N atoms. Massive 
MS stars is supported against gravity with hydrogen burning 
via CN-cycle reactions. With lower CN abundance, higher 
temperature is required to supply enough to support the star, 
thus leading to smaller MS radius. 

The lower panel of Figure[3]shows the protostellar evolution 
with accretion rate ten times higher than the fiducial value. 
The mass-radius relation varies with metallicity in the same 
manner as in the fiducial cases. Compared to the fiducial case 
at the same metallicity, a protostar has larger radius and ar- 
rives to the ZAMS at higher stellar mass, which is consistent 
with the dependence on accretion rates discussed above. Clear 
differences are found only for Z = and 10"^ Z©. The proto- 
star turns to inflate during the KH contraction in these cases. 
We separately focus on these cases in § IB.l.ll below. 

3. LIMITS OF MASS ACCRETION 
3.1. Radiative Feedback Effects 



Protostellar luminosity increases with growth of the star by 
accretion. Radiative feedback effects will eventually termi- 
nate the accretion. Below we examine the effects of radiation 
pressure exerted on a radiative precursor and dust cocoon, 
as well as of the expansion of an H II region. Our calcula- 
tions have already included radiation pressure on the radiative 
precursor We examine the other effects by separate analytic 
treatments and evaluate when these effects limit the accretion 
as a function of metallicity. 

3.1.1. Radiation Pressure on Radiative Precursor 

A protostar growing with very high accretion rate exerts 
strong radiation pressure on a radiative precursor, which is 
a optically thick part of the gas envelope. This effect is 
found in the calculated protostellar evolution with lOMfid for 
Z < 10"^ Zq (the lower panel of Fig. |3]l. In these cases, 
at ~ 60 Mq the KH contraction turns to the abrupt ex- 
pansion (also see Omukai & Palla 2003). The expansion oc- 
curs because radiation pressure decelerates the accreting flow 
before reaching the stellar surface. The stellar surface pres- 
sure, which balances the ram pressure by the flow, becomes 
too low to maintain the hydrostatic structure of the star. The 
expansion occurs when the total stellar luminosity Ltot ap- 
proaches the Eddington limit LEdd = 47rcGM*/Ke, where 
is the Thomson scattering opacity. Although the luminosity 
from stellar interior remains below the Eddington limit, 
the total luminosity Ltot from the protostar has an extra contri- 
bution from the accretion shock: 



Ltot ~ ^* "t" ^ai 



GMM 



R^, 



(7) 



Note that the internal to accretion luminosity ratio L^/L^^c is 
equal to the timescale ratio faee/fKH- Since Lace (xM/R^:, Lace 
becomes very high and Ltot reaches the Eddington limit for a 
high enough value of M during the KH contraction, where the 
radius approaches the ZAMS value. This happens for the ac- 
cretion rate exceeding the critical value M^ ~ 4 x lO"-' Mq/yy 
for the case of Z = and depending on metallicity only weakly 
(Omukai & Palla 2003; Hosokawa & Omukai 2009). This is 
indeed the case for Z < 10"^ Z© with lOMgd. In such cases, 
the protostars cannot keep accreting at its original rate after 
the abrupt expansion. If the accretion terminates at this mo- 
ment, the final stellar mass will be ~ 80 M© (Fig. |5]l. 

3.1.2. Radiation Pressure on Dust Cocoon 

We here estimate the effect of radiation pressure on the 
dust cocoon. First, we consider how the radiation propa- 
gates outward in the accretion envelope. Most of the stellar 
radiation is initially emitted in UV and optical wavelengths. 
The stellar radiation travels freely between the photosphere 
and dust destruction front. Subsequent outward propagation 
in the dust cocoon depends on whether the cocoon is opti- 
cally thick or thin for the stellar UV/optical light. If the co- 
coon is optically thick, most of the direct light is absorbed 
near the dust destruction front and is re-emitted in infrared 
wavelengths, then propagates outward diffusively. This is the 
case for the present-day protostars still embedded in their na- 
tal dense cores, which are observed as bright infrared sources. 
With lower abundance of dust grains, however, the cocoon be- 
comes transparent even for the stellar UV/optical light. In this 
case, the direct stellar light escapes from the cocoon without 
absorption and re-emission. Figure|2]schematically shows the 



Protostellar Evolution of Low-metallicity Stars 



5 




Fig. 4. — The stellar upper mass limit by radiative feedback at the fiducial accretion rate as a function of metallicity. The thin dotted and dashed lines shows 
two characteristic epochs in protostellar evolution; the maximum of the stellar radius (labeled as max') and arrival to the ZAMS('ZAMS'). The dust cocoon 
is optically thick in the upper-left hatched area and thin in the rest. The right gray-shaded areas denote the region where the mass accretion is prohibited either by 
the radiation pressure on the dust cocoon ('Prad '): the expansion of an H II region ('HII region'), or the stellar lifetime ('lifetime'). The characteristic mass-scales 
of prestellar cores by fragmentation are also presented. The thick-dashed and dot-solid lines represent the cores produced as a result of fragmentation induced by 
dust cooling and line cooling, respectively. 




Fig. 5. — Same as Fig|4]but for the accretion rate ten times higher than the fiducial value. The right-bottom gray-shaded rectangular area denotes the region 
where mass accretion is prohibited by the radiation pressure exerted on the radiative precursor 
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Fig. 6. — Evolution of the dust cocoon at different metallicities; Z = 
Zq (solid lines), 10"^ Zq (dashed lines), and 10^ Zq (dot-dashed lines). 
The fiducial accretion rates are adopted at each metallicity. Upper panel: 
Position of the dust destruction front, i.e., the inner edge of the dust cocoon. 
Lower panel: Optical depth of the cocoon for the stellar UV and optical light 
Tc. The cocoon is optically thick if Tc > 1, which is indicated by the dotted 
background. 

Structure of the flow for these two cases. The optical depth of 
the cocoon is roughly estimated as 



'«UV/optPd^d, 



(8) 



where Kuv-opt is dust opacity for optical and UV light, and the 
suffix "d" denotes physical quantities at the dust destruction 
front. We regard the cocoon as optically thin (thick) if Tc < 1 
(> Irrespectively). The optical depth can be evaluated as 
follows. The position of the dust destruction front is 



Ri ~ 35 AU 



105 Lr. 



1/2 



(9) 



(Wolfire & Cassinelli 1987). The upper panel of Figure |6] 
shows evolution of the dust destruction front at Z = 10"**, 10"^, 
and Zq for the fiducial accretion rates. The front moves out- 
ward with increasing stellar mass as the total luminosity Ltot 
increases. We adopt a typical dust opacity in the optical range 
(e.g., McKee & Ostriker 2007), 



^uv-opt ■ 



:200cm7g( — ). 



For a free-falling flow, the density at is 



Pd = 



M 



(10) 



(11) 



A-nJlGM^Rl 



We derive Tc c>c pd^d oc L J/"^ by substituting equations (|9]l and 
(fTTT i into (O. Hence the optical depth decreases with the in- 
crease of stellar mass and luminosity. Even if the cocoon is 
initially optically thick, it becomes optically thin at some mo- 
ment in the evolution. Evolution of the optical depth shown 
in lower panel of Figure |6] demonstrates this transition. At 
Z = 10~^ Zq, for example, the cocoon becomes transparent 
for the direct fight at ~ 3 Mq. Note that the optical depth 
Tc cx Kuv/optPd ZM does not linearly scale with Z, owing 
to another dependence on M. Figures |4] and |5] summarize the 
evolution of the dust cocoon for all cases. 

Next, we consider the condition for the flow to overcome 
the radiation pressure for cases where the cocoon is optically 
thick or thin, respectively. In the optically thick case, most of 
the direct radiation is absorbed in a thin layer near the dust 
destruction front and isotropically re-emitted subsequently in 
infrared wavelengths. Hence the absorbing layer receives out- 



ward momentum flux of 



For the accretion to proceed. 



the inward momentum flux of the flow must exceed this out- 
ward flux (Wolfire & Cassinefii 1987), 



PdMd > 



AttRIc 



where is the free-fall infall velocity. 



2GM, 
Ra 



(12) 



(13) 



Using the continuity equation M = AnR^pu, equation (fT2] i 
leads to a condition for accretion rates. 



M> 



(14) 



Additionally, the flow in the cocoon must overcome the radi- 
ation pressure by the infrared re-emission. This condition can 
be stated as that the stellar total luminosity be less than the 
Eddington limit defined with dust opacity for infrared light 
Kir: 

_AttcGM, 

i-tot < i-Ed.IR = ■ (Ij) 

Kir 

Otherwise, radiation pressure decelerates the accretion flow 
before reaching the dust destruction front. The infrared opac- 
ity in the local star-forming regions takes its maximum at 
Tdust =i 600 K (e.g., Poflack et al. 1994). We adopt the Rosse- 
land mean opacity at this maximum as the infrared opacity 



Kir ~ 8 cmVg 



(16) 



On the other hand, for an optically thin cocoon, the stellar 
radiation is not converted to infrared radiation by absorption 
and re-emission. The only condition is that the total stellar 
luminosity be less than the Eddington hmit for the opacity at 
the UV/optical light Kuvopt. 



Ltot < Lfid, 



47rcGM, 



uv-opt ■ 



(17) 



^uv-opt 

Figure s|4] and |5] show the limit when the radiation pressure 
becomes so strong to violate the above conditions. For Z = Zq 
with the fiducial accretion rate, for example, the flow is opti- 
cally thick until M^, ~ 30 Mq. For an optically thick flow, 
equation (fT4l l usually put the stricter condition than equation 
(fISl l. and the upper limit falls on ~ 8 Mq for Z = Zq (Fig- 
ure IH. At Z = 10"^ Zq, the optically thick/thin transition of 
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the flow occurs at ~ 3Mq and the accretion continues until 
M, > 50 Mq, where the condition ([TtI i is violated. Both fig- 
ures show that the limiting mass increases with decreasing 
the metallicity. Radiation pressure on the dust cocoon does 
not limit mass accretion at Z < lO"-' Zq. 

3.1.3. Expansion of an H II Region 

The stellar UV luminosity increases with its mass. Strong 
UV photons ionize gas near the star and an H II region de- 
velops inside the accreting envelope. For the free-falling flow 
with accretion rate M, the radius of the H II region Run can be 
expressed as a function of the UV photon number luminosity 
S: 

RHii = RphSxp(J-^, (18) 

where /?ph is radius of the photosphere and the critical UV 
luminosity Scr is defined as. 



Srr — 



aM 



(19) 



where a is the recombination efficiency and /iH is mass per 
hydrogen atom. As long as 5 < S^, the HII region is con- 
fined near the stellar surface. Once the stellar UV luminosity 
exceeds 5cr, the HII region grows exponentially and soon en- 
closes all materials in the natal core. Figure |7] presents the 
evolution of the protostar and its H II region for Z = 10"'* Zq 
with the fiducial accretion rate. In this case, the photosphere 
locates beyond the accretion shock front at the stellar surface 
except short duration around 2± 6M^. The UV photon 
number luminosity is calculated by 



hv 



(20) 



where t'L is the Lyman-limit frequency, h is Planck constant, 
and B{v, Teft) is Planck function with the effective temperature 
Teff. Figure [7] shows that the H II region begins to expand 
exponentially at M, ~ 40 Mq. 

The evolution of the HII region and the accretion flow ex- 
periences the following two stages, depending on the ratio of 
the radius of the HII region /?hii and the gravitational radius 
of an ionized gas, 

R^ = ^^^ (21) 

<^HU 

where chii is the sound speed of ionized gas and 0Ed = 
1 -Ltot/i'Ed is a correction factor taking account of repulsive 
force by radiation pressure. When Run < Rg, the accreting 
flow supersonically passes through the ionization front. Al- 
though the temperature and pressure jump up behind the ion- 
ization front, this hardly influences the flow dynamics. The 
heated gas inside the H II region is also gravitationally bound 
to the central protostar. Once Run exceeds Rg, the accreting 
flow becomes subsonic and dynamical expansion of the H II 
region disrupts the accreting envelope. The gas inside the H II 
region begins to move outward owing to the pressure excess 
of the ionized gas. A shock front emerges ahead the ioniza- 
tion front, and sweep up materials into a shell around the H II 
region. 

Strictly speaking, the density distribution of the flow de- 
viate from the free-fall one even without rotation. Radia- 
tion pressure via photoionization decelerates the flow and in- 
creases the density. Omukai & Inutsuka (2002) showed that 
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Fig. 7. — Evolution of the protostar and surrounding H II region at the 
fiducial accretion rate for metallicity 10^ Zq. The thick and thin solid lines 
represent the position of the stellar surface, i.e., the accretion shock front, 
and the photosphere within the accreting flow. Positions of the ionization 
front Shii and gravitational radius for an ionized gas Rg are presented with 
dashed and dotted lines. The steady mass accretion is disturbed by dynamical 
expansion of the H II region when Rhii > ^g, which is denoted by the gray 
shaded background. The dot-dashed fine represents the dust destruction front, 
the region above which coiresponds the dusty accreting envelope. 

this effect can significantly delay the dynamical expansion of 
an H II region around the first star (Z = 0). However, McKee 
& Tan (2008) argued that relaxing the spherical symmetry re- 
duces this effect. In the rotating core collapse, density in the 
polar region decreases by a large factor by the epoch when 
an H II region forms. As a result, radius of the H II region 
reaches Rg before the trapping effect becomes important (see 
§ 5.2 of Mckee & Tan 2008 for detail). The trapping effect be- 
comes even weaker for lower accretion rates, i.e., for higher 
metallicity. In this paper, although we treat the spherically 
symmetric flow in estimating the radiation feedback, small 
asphericity is assumed to exist to alleviate the trapping effect 
by the photoionization. 

Figure |6] shows that the dynamical expansion occurs at 

~ 100 Mq at Z = 10"^ Zq. The steflar mass limit by the 
HII region expansion is summarized in Figures |4] and |5] as a 
function of metallicity. With the fiducial accretion rates, the 
expansion of an H II region is the primary feedback effect 
disrupting the accreting envelope for Z < 10"^ Zq. The mass 
limit increases from ~ 30 Mq to 300 Mq with decreasing 
metallicity. This is because higher accretion rate can quench 
the H II region of the more massive star. With ten-times the 
fiducial accretion rate, the mass limit exceeds 1000 Mq for 
Z < lO"'* Zq. In these cases, not the expansion of an HII 
region, but the radiation pressure on th e radia tive precursor 
limits mass accretion for Z < 10"^ Zq (§ 13.1.11 ). 

The mass limit by the H II region expansion is ~ 
300 Mq for Z = with the fiducial accretion rates. This 
happens to be comparable to that by McKee & Tan (2008), 

~ 100-300 Mq. However, they consider mass accre- 
tion via an accretion disk formed by rotating core collapse. 
The breakout of an H II region beyond Rg occurs earlier at 

~ 50-100 Mq. Photoevaporation of the accretion disk 
finally quenches mass accretion. The same process can work 
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in halting the accretion also for finite metaUicities, and should 
be explored in future work. 

3.2. Mass Limit by Stellar Lifetime 

We have shown that mass accretion to a protostar is limited 
by radiative feedback effects, i.e., radiation pressure exerted 
on a dust cocoon for Z > 10"^ Z©, and expansion of an H II 
region for lower metallicity. Figures |4] and |5] show that these 
limits are reached after the protostar arrives at the ZAMS. For 
very high feedback limits, accreting stars might end their lives 
before reaching such limits. The maximum mass of the star 
can gain in its lifetime is roughly given by 

~ / Mdt (22) 
Jo 

where fMS is the lifetime of a ZAMS star of M,. This limit is 
also presented in Figures |4] and |5] Here, we adopt Tms of pri- 
mordial stars from Marigo et al. (2001) forM, < lOOM© and 
from Schaerer (2002) forM, > 100 M©, as their dependence 
on metallicity is weak. We see that the limit by the lifetime 
is less stringent than those by the radiative feedback effects. 
Hence the final stellar mass will be set by the feedback effects, 
as long as the prestellar core is massive enough. 

4. CHARACTERISTIC CORE MASS SCALES 

The mass limit derive above is reached only in the case that 
the mass reservoir, i.e., the mass of the natal core, is sufficient. 
Here, we discuss the characteristic mass of the prestellar cores 
and compare them with the mass limits by the feedback. The 
cores are formed as a result of fragmentation of a more mas- 
sive cloud. According to analytical (Schneider et al. 2002; 
Larson 2005) and numerical (Bromm, Coppi & Larson 2002; 
Jappsen et al. 2005; Clark, Glover & Klessen 2008; Smith 
et al. 2009) arguments, fragmentations likely occur around 
temperature minima of thermal evolution tracks (see Figure 
[U, and cores of about the Jeans mass Mj are produced. This 
figure shows that the tracks at 10"^ Zq < Z < 0.1 Zq expe- 
riences two local temperature minima, i.e., one at lower den- 
sity 10^"^cm"\ where Mj ~ 100- 1000 Mq, and the other 
at higher density 10^"''^cm"^, where Mj ~ 0.01-0.1 Mq. 
The former is induced by the line cooling, i.e, cooling by 
the H2 lines for < lO'^Z©, and the metal fine-structure lines 
for higher metallicity, while the latter is by the dust cooling. 
In clouds with extremely low-metallicity, dust-induced min- 
imum disappears, thus only high-mass fragmentation is al- 
lowed. On the other hand, at Z ~ Z©, two minima almost 
merges and result in a single characteristic mass-scale. The 
typical core mass in the cases with two fragmentation epochs 
is not clear yet. We speculate that massive and small cores 
corresponding to two fragmentation epochs exist simultane- 
ously, and their number ratio depends on the efficiency of 
dust-induced fragmentation. 

Those characteristic core mass-scales are also shown in Fig- 
ure |4] and |5] for comparison with the stellar mass limits by the 
feedback. The characteristic core masses by the dust-induced 
fragmentation is much smaller than the stellar mass limits. In 
such cores, the final stellar masses are bounded by the amount 
of mass reservoir rather than the stellar feedback. On the 
other hand, the typical core masses by the line-induced frag- 
mentation are comparable to the mass limit by the radiative 



effects. In such massive prestellar cores, therefore, the final 
stellar masses should be limited by the radiative effects and 
very massive stars can be formed in low-metallicity environ- 
ments. 

5. SUMMARY 

We studied protostellar evolution as a function of the metal- 
licity, using the results to evaluate the significance of radiative 
feedback effects. For metaUicities between zero and Solar, we 
showed how the feedback terminates the mass accretion and 
determines the final stellar mass. Our results are summarized 
as follows: 

First, we estimated accretion rates onto a protostar from 
the thermal evolution of a prestellar core at each metallicity. 
The accretion rates are higher at the lower metallicity, ranging 
from M ~ 10"^ Mg/yr at Z = to ~ 10"^ Mg/yr at Z = Zq. 
This is because the core temperature is higher at the lower 
metallicity. 

Next, we numerically calculated the protostellar evolution 
with the derived accretion rates at each metallicity. The calcu- 
lated evolution depends on the metallicity. For lower metal- 
licity, a protostar has larger radius and arrives at the ZAMS at 
higher stellar mass. These variations reflect the difference in 
the accretion rates. 

We also examined radiative feedback effects that may halt 
the mass accretion onto a protostar, i.e., radiation pressure on 
a radiative precursor and dust cocoon, and expansion of an 
H II region. The radiation pressure exerted on the dust co- 
coon limits mass accretion for Z > lO""* Zq. The upper mass 
limit increases with decreasing metallicity. At low metallic- 
ity, the high accretion rate and low dust opacity enable the 
accreting flow to overcome the radiation pressure barrier. For 
Z < 10"^ Zq, the expansion of an H II region instead will dis- 
rupt the accretion envelope. The mass limit resulting from the 
H II region expansion also increases with decreasing metallic- 
ity because higher density in the flow hinders the development 
of the HII region. These feedback effects become significant 
after the protostar reaches the ZAMS, but before the stellar 
lifetime is over. With the accretion rates ten times higher 
than the fiducial values, radiation pressure on the radiative 
precursor limits the steady mass accretion at ~ 100 Mq 
forZ< 10"^ Zq. 

Cores of two typical mass-scales are produced by fragmen- 
tation over wide ranges in metallicity. One scale comes from 
line cooling (100- IOOOMq), and the other from dust cool- 
ing (0.01 - IMq). The mass limit due to the radiative feed- 
back effects is far higher than the typical masses with dust 
cooling, and comparable to those with line cooling. In mas- 
sive cores formed by line-induced fragmentation, the radia- 
tive feedback effects terminate mass accretion and set the final 
stellar masses. 
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